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I. INTRODUCTION 

In its original form, classical dynamical density func- 
tional theory (DDFT) was derived by Marconi and Tara- 
zona [T] in 1999 for spherical, i.e., isotropic, colloidal 
particles. Their derivation started from the Langevin 
equation for spherical particles that interact via a 
pair potential. Later, in 2004, DDFT was rederived by 
Archer and Evans (3] from the Smoluchowski equation 
that corresponds to the Langevin equation for interact- 
ing spherical particles. In 2007, DDFT was generalized 
by Rex, Wensink, and Lowen [3] to systems of uniaxial 
anisotropic particles with orientational degrees of free- 
dom. This generalization is based on the Smoluchowski 
equation for rigid rods [5|. It made DDFT applicable to 
the important class of uniaxial liquid crystals. 

Nowadays, it is already possible to produce colloidal 
particles with rather complicated shapes including biax- 
ial particles. Although static classical density functional 
theory (DFT) has presently available very powerful tools 
like fundamental measure theory [B] that allow to con- 
sider also such complicated colloidal particles in the con- 
text of DFT, the dynamics of these biaxial particles could 
up to now not be investigated on the basis of DDFT. For 
these reasons, it is of high importance to push forward 
the development of DDFT. 

In this paper, we present a further generalization of 
DDFT, which is now also applicable to biaxial particles. 
This extension of DDFT contains the previous DDFT 
equations as special cases and does not assume a cer- 
tain shape for the colloidal particles. Instead, it is de- 
rived for arbitrarily shaped colloids. In comparison with 
the former DDFT approach, this leads to three indepen- 
dent rotational diffusion coefhcients instead of only one. 
Since our new DDFT equation holds also for screw-like 
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particles, it takes even a possible translational-rotational 
coupling into account. Additionally, we consider a possi- 
ble self-propulsion mechanism of the particles so that our 
results are also relevant for the investigation of the collec- 
tive dynamics of active particles like swarms of swimming 
microorganisms as, for example, protozoa [7]. 

The paper is organized as follows: after giving a short 
overview in Sec. |ll] about anisotropic colloidal particle 
shapes that can already be synthesized, we present our 
derivation of the extended DDFT equation in Sec. |III| 
and discuss special cases that are known from literature. 
Sec.|lV]is addressed to possible apphcations of the DDFT 
equation. Finally, we give conclusions and mention pos- 
sible further extensions of DDFT in Sec. |Vl 



II. GEOMETRIC CLASSIFICATION OF 
COLLOIDAL PARTICLES 

Induced by technological advance in the processing 
of nanomaterials, a large number of differently shaped 
colloidal particles became synthetizable during the last 
years. The different shapes of these colloidal particles 
can be classified by means of their geometric properties. 
Figure [T] shows a detailed classification of colloidal shapes 
with respect to symmetry and convexity. Such a classifi- 
cation is of big importance since colloidal particles may 
form a huge set of mesotropic phases (mesophases) [H |9] 
that go along with different states of translational and 
orientational order. The possible states of translational 
and orientational order depend strongly on the shapes of 
the particles and a classification of their shapes is there- 
fore also a classification of the possible phases that these 
particles may evolve. 

The most simple and at once full symmetric, i.e., 
isotropic, shape is the sphere. This is the traditional 
shape for colloids in theoretical soft matter physics, be- 
cause it is simple to produce and due to a lack of orien- 
tational degrees of freedom relatively simple to describe 
theoretically. Since spheres possess only translational de- 
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FIG. 1. Classification of synthetizable colloidal particles with respect to their shape. Geometrical properties that were used to 
classify the shapes are symmetry and convexity. 



grees of freedom, they solely appear in the completely dis- 
ordered isotropic phase and in the crystalline state [TU] . 
The shape of a sphere is globally convex and there is no 
non-convex analog with full symmetry. All other colloidal 
particle shapes are anisotropic and either uniaxial or bi- 
axial. The characteristic property of uniaxial particles 
is a symmetry axis, whose orientation is denoted by the 
unit vector u in the following. These particles have rota- 
tional symmetry and possess one orientational degree of 



freedom in two spatial dimensions and two orientational 
degrees of freedom in three spatial dimensions. Uniaxial 
particles are further distinguished into apolar and polar 
particles. An uniaxial particle is called apolar, if it has 
head-tail symmetry and polar otherwise. Rod-like parti- 
cles [TTJ[T2] like spherocylinders, spheroids, and ellipsoids 
are the most simple anisotropic colloidal particles. They 
are convex and apolar and of big importance since they 
may evolve the industrially important nematic phase and 
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serve as excellent model systems for most liquid crystals 
[13H15j . A further member of convex and apolar particles 
are the platelets [TSHSU]. They have a similar phase dia- 
gram to rod-like particles with a strong affinity to form 
columnar stacks ^21j . Systems of such disk-like particles 
are realized in nature for example by clay suspensions 
Examples for non-convex apolar particles are 
dumbbells (dimers) [25 -27 , that are produced by mer- 
gence of two spheres of equal size, and rings [23 [H], that 
can be made by etching from colloidal spheres that are 
partially embedded in a metal layer. The complement 
of apolar particles is built by the polar particles, that 
have no head-tail symmetry. A famous member of this 
particle class are the Janus particles [30H32] . They are 
spheres with a different coating at one half of the sur- 
face. The original Janus particles had a hydrophilic and 
a hydrophobic coating. Nowadays, one coating is often 
reactive like a platinic coating that decomposes hydrogen 
peroxide catalytically. Such particles are immersed into 
a hydrogen peroxide solution to realize active particles 
(micro-swimmers) that are driven by an intrinsic drive 
[33j . Cones are a further member of uniaxial polar par- 
ticles. Carbon nanocones appear naturally in graphite 
and do not need to be produced by an elaborate 
method. By the mergence of two spheres with different 
diameters, one obtains a pear-like particle [371 135] ■ Pears 
and also bowls [521 SO] are non-convex particles that are 
uniaxial and polar. The latter stack into each other and 
form columnar structures [41) . 

Particles with less symmetry are biaxial. They are the 
complement to the uniaxial particles in the class of the 
anisotropic particles. Biaxial particles have either only 
discrete symmetries, like inflection symmetry and dis- 
crete rotational symmetry, or are completely asymmet- 
ric. In both cases, the biaxial particles have three orien- 
tational degrees of freedom and a unit vector is no longer 
sufficient to describe their orientation. Instead, two per- 
pendicular unit vectors or Eulerian angles have to be used 
[12]. Due to the additional orientational degree of free- 
dom, the phase diagrams of biaxial colloidal particles are 
much richer than those for uniaxial particles |43j . Convex 
colloidal particles with discrete rotational or inflection 
symmetry are, for example, polyhedra like cubes [I^HIT] 
and tetrahedra 03 [49], boards [50], pyramids [5TH53] . 
and regular patchy particles [5iH57] . The latter differ 
from Janus particles by a patchy coating with a regular, e. 
g., tetrahedral, arrangement. Non-convex particles with 
discrete rotational or inflection symmetry include special 
colloidal molecules that are realized by more than two 
spheres that are merged in a regular arrangement. Exam- 
ples for this include triniers [SS] consisting of three equal 
spheres and chiral particles [59', (60] consisting of many 
equal spheres in a helical arrangement, multipod-shaped 
nanocrystals [HIl 112 , stars [351 113] ) and some lock-and- 
key particles [64] . Patchy particles may also belong to the 
class of colloidal particles without any kind of symmetry. 
This is the case, if the patches are arranged or sized in an 
irregular way. Irregular patchy particles that are made 



by coating of spherical particles are always convex. Col- 
loidal molecules of arbitrary shape and size belong on 
the other hand to the completely asymmetric colloidal 
particles that arc not convex [55HB5] . 

III. DERIVATION OF THE DDFT EQUATION 

In this derivation, we consider a set of N asymmetric 
rigid particles in a solvent with dynamic (shear) viscos- 
ity 77 and neglect possible additional (for example vibra- 
tional) degrees of freedom. We choose the center-of-mass 
positions fi = {xi^i,X2,i,X3,i) and the Eulerian angles 
mi — (ipi, 0i, Xi) with i = 1, . . . , A'' to describe their po- 
sitions and orientations completely. Alternatively, the 
orientation of the particles could also be described by 
means of two perpendicular axes |69| . but for our pur- 
poses, the use of Eulerian angles is more appropriate, 
since they do not involve additional geometric constraints 
and lead to simpler equations with a more compact nota- 
tion. The angular velocities cJ; that describe the instanta- 
neous rotational motion of the particles can be expressed 
in terms of the Eulerian angles and their temporal deriva- 
tives |70| . For convenience, we use the convention of Gray 
and Gubbins [7T], which is equivalent to the second con- 
vention of Schutte [70], for the Eulerian angles, since with 
this convention, the first two Eulerian angles (f> and 6 
are identical to the usual azimuthal and polar angles of 
the spherical coordinate system, respectively. The whole 
set of particles is then characterized by the positional 
and orientational "multivcctors" = (ri,...,r/v) and 
— {wi, . . . ,wn), respectively. For completeness, we 
also introduce the abbreviation uj^ — (wi, . . . , ujn), here. 
The particles are exposed to the (time-dependent) total 
potential 

C/(r^, t) = U,^t ir", t) + (7i„t (r^, t5^) , (1) 
which consists of the external potential 

N 

C/e,t (r ^, t^"^, = 5Z (2) 

i=l 

and the total particle interaction potential 

N 

Uint{r^,-^^) = ^ U2{ri,rj,Wi,Wj) . (3) 
i<j 

For both the one-particle interaction potentials 
Ui{fi,-Sji,t) and the two-particle interaction poten- 
tials U2{ri,rj,Wi,Wj), we assume pairwise additivity. 
Moreover, we neglect many-particle interaction poten- 
tials of higher order than pair interaction potentials. We 
further introduce the A'^-particle probability distribution 
function P(r^, w^, t) for the probability density to find 
the N particles at time t with the orientations at 
the positions r^. Successive integration of this function 
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with respect to its positional and orientational degrees 
of freedom leads to the n-particle density |3J 



(4) 



X /dr!„+i---/dr!ArP(r^,t^^,t) 
Js Js 



where V = and 5 = [0, 27r) x [0, tt) x [0, 27r) are the 
domains for spatial and orientational integration, respec- 
tively, dV = dxidx2da;3 and dfl — d(j)d9 sm{9)dx are the 
corresponding differentials, and 
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dV = / dxi I dx2 I dx3 , 

V JQ Jo Jo 

dn= d4> d6isin(6') dx 
; Jq Jo Jo 



(5) 



are common abbreviations. 



A. Smoluchowski equation 

We start with the derivation of the Smoluchowski equa- 
tion for the overdamped Brownian dynamics of N self- 
propelled biaxial particles. In analogy to the uniaxial 
passive case (see reference jS]), we define the translational 
gradient operator = (dx^ , , dx^) and the rotational 
gradient operator — ih, where i is the imaginary unit 
and L = (Lxi-,^x2t^x3) is the angular momentum oper- 
ator, which can be expressed in terms of the Eulerian 
angles by [7D] 



iLx 



d d 
■ cos(0) cot(6') — - sm(0) — 



■ cos((/)) csc(6') 



d_ 

dx ' 



d d 
■s\n{4>) cot{0)—+cos{(j)) — 



(6) 



-I- sin(0) csc(0) 
d_ 



dx 



We further define the vectors = {f^,w^) and 
(r^, w^) with — df^/dt and the operators V^n 
(V,-=-i,...,Vf„), V^jv = (V^i, • . • , Vs„), and Vjjv 
(VpN, V^jv) and write down the continuity equation 



d_ 
di 



P(y^,<) = -V,-«.(o^P(y~,t)), 



(7) 



which is a trivial generalization of the continuity equation 
for passive rods that is described by Dhont in Ref. ^ . On 
the Brownian time scale, the total force and torque, act- 
ing on an arbitrary particle i G {1, . . . , iV} are zero. The 
total force and torque consist of the force f'i^^ (y^, t) and 



torque T^^^ (y^, t) due to the activity of the self-propelled 
particle i, the hydrodynamic force fI^\i^) and torque 



t'>^\i^), the interaction force F^'''{f^,t) and torque 
fP{f'^,t) due to the potential U{f^,t), and the Brow- 
nian force F^^''\f'^ , t) and torque Tj;^'^\f'^,t). With the 
definition X = {Xi,...,Xn) for X e {F^ ' f ( ' ), i?( ' )} 
and the abbreviations 



i?m(y^,t) = (FW(r^,t),fm(y^,i)), 



(8) 



this force balance for the N colloidal particles can be 
expressed by 



= i?(^)(y^,t)+i?(H)(y^) 
-l-i?m(y^,t) + i?(i^'-)(y^,t) 



(9) 



The forces and torques resulting from the self-propulsion 
mechanism of the particles are supposed to be constant 
with respect to their orientations in the respective body- 
fixed coordinate systems, but their strengths may vary 
slowly with time. We denote these forces and torques 
for a certain particle i in body-fixed Cartesian coordi- 
nates by the vector i?o^^ iri,t) for i = 1, . . . , iV and the 
corresponding vector in space-fixed coordinates by 

Kl^'>{n,w,,t) = no\m,)0of{n,t) (10) 
with the diagonal block rotation matrix 

TZq{w) = diag(R(ro),R(t5)) , (11) 
where the rotation matrix R(to) is defined by 



(12) 



(13) 



R(t5)-R3(X)R2(^?)R3(0), 

R-i(t5) =RT(t5) =R3(-0) R2(-e)R3(-x) 
with the elementary rotation matrices 

(cos{(f) — sin((/3)^ 
1 
sin((p) cos(<y9) 

(cos{tp) sm{tp) 0^ 
-sin(93) cos{ip) 
1, 

^ (A) 

Note that Kq / depends most often only on time t, but 
one could also think of swimming microorganisms in a 
poisoned environment, where ^q^^ depends on r^. 
To simplify the notation in the following, we collect all 
the N vectors K[^\ri,Liji,t) in the vector 



i?(A)(y^,i) = 7^-l(t5^)i?(^^^(f^,^) 



(14) 
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with the 6N x 6 A'^-dimensional rotation matrix 

n{m^) = diag(7^o(t5l), . . . , noiwN)) (15) 

and the 6iV-dimensional vector 

^ {0,'^{n,t),...,K^,'^ir^,t)) . (16) 

Next, we focus on the hydrodynamic force and torque. 
They are given by 



(17) 



with the microscopic friction matrix [S] 

"^(y^) — ( tKT^n} -r-BR^Ni ] ; (18) 



where T^^if"), T™(f^), T'^'^if^), and T^^if^) are 
3N X 3iV-dimensional submatrices. The submatrices 
YTTj-yW-j g^^^ YRR^yW-j correspoud to pure translational 
and rotational motion, respectively, while xTR^yW-j 
j-KTfjN-^ have to be taken into account for particles 
with a translational-rotational coupling as, for example, 
screw-like particles. For many other particles like those 
that are orthotropic, however, T™{f^) and T^'^(y^) 
vanish. In the following, we neglect hydrodynamic in- 
teractions between the colloidal particles. With this as- 
sumption, the microscopic friction submatrices simplify 
to the block diagonal matrices 

TTT(^a') = diag(TTT(^^), . . . , T^T (^^)) , (19) 

T™(^^) = diag(T™(mi), . . . , T™ (^^,)) , (20) 

T«T(^A') = diag(Tf;^(^,), . . . , TYNi^^)) , (21) 

r^^im^) = diag(Ti^i^(^i), . . . , TTNi^N)) (22) 

with the 3 x 3-dimensional matrices 

T^im,) = 77R-i(^,)KR(^5.) , (23) 

T™(w,) = 77R-i(t5,)C(')^R(t5.) , (24) 

T«T(m,) = r^R-'{m,) C^^) R{w,) , (25) 

- r,R-i(t5,) R(t5,) (26) 

for i = 1, . . . , iV, which are related to the translation ten- 
sor K, the coupling tensor C*-^-*, its transposed C^^^""", and 
the rotation tensor fl^^^ [72] by an orthogonal transfor- 
mation with the rotation matrix R(ro). The tensors K, 
C(S), and 

are constant and depend on shape and 
size of the colloidal particles that are considered, but are 
independent of the viscosity of the solvent. In addition, 
C'^^ and r2*^^^ depend also on the reference point S, for 
which the center-of-mass position of the considered col- 
loidal particle should be chosen. In the special case of no 
hydrodynamic interaction, the inverse of the microscopic 
friction matrix 



/3V{f 



(27) 



with the inverse temperature /? = 1/(A:bT), the Boltz- 
mann constant fee , and the microscopic short-time diffu- 
sion matrix 



1 pRT(y.iV) 2?RR(y-A^) 



(28) 



which we need in the following, has the same structure 
as the microscopic friction matrix. We further have the 
equation 



(29) 



for the interaction force and torque. Moreover, the Brow- 
nian force and torque F'^^'\f^,t) and f^^'\f^,t) can be 
derived from the equilibrium condition 



lim t) cx e 



(30) 



when (y^, t) is neglected and the vector in Eq. 
([t]) is expressed in terms of the vectors , K'--^^ (y^, t) , 
and i?(^')(y^,t) with the help of Eq. This results in 



(31) 



Using Eqs. (|9|, ([T7|, ([29]), and (|3l]), the Smoluchowski 
equation 



d_ 
Ft 



Pir,t) = cpir,t) 



(32) 



(33) 



with the Smoluchowski operator 

£-V^«-(l?(y^)-(/3Vj-«C/(f^,t) 
-/3i?(A)(y^,t)-f Vj-«) 
follows now directly from the continuity equation Q 

B. DDFT equation 



Next, we proceed in our derivation by applying the 
integration operator N J^V2- ■ ■ J^Vn J^^2- ■ ■ J^^n 
from the left on the Smoluchowski equation ([32]) and ob- 
tain the expression 



p{h t) = Vj- • (2?(y) • ( - PKaH t)p{l t) 



d_ 

+V^p{l t) + Pp{l t)V^U,{l t) - pK{l t)j) 



with the short-time diffusion tensor |73| 



(35) 



for the one-particle density p{f,t) = p^^-'(y, <), where we 
omitted the index 1 in fi and Wi and used the abbrevi- 
ations y = (f,73), = (V^., Vs), KAiht) = K[^\ht), 



6 



and K(j,t) = {F{f,t),T{f,t)). When we further intro- 
duce the integration operator 



d93 = dv an 

Jv Js 



(36) 



with the total integration domain & — V x S and the 
corresponding differential d5J = dV^dfJ, the average force 
F{f, t) and torque T(y, t) due to the interaction with 
other particles in Eq. ( 34 1 are given by 



Kilt) 



(37) 



In equilibrium with Kp^{f,t) — and Ui — Uiij), Eq. 
( [M] ) reduces to the first equation of the Bogoliubov- 
Born-Green-Kirkwood-Yvon hierarchy for molecular flu- 
ids [7T]: 



(38) 



Here, a zero in the index of a function denotes the 
time-independent equilibrium state of this function. For 
example, the function po(yl denotes the equilibrium 
one-particle density field that corresponds to the time- 
independent prescribed external potential Uiij). On the 
other hand, we have in equilibrium the relation 



VfPo(r)+/3po(r)Vf[/i(y) 



'^■^exc[Po(y)] 

<5po(y) 



(39) 



with the equilibrium Helmholtz excess free-energy func- 
tional ^cxc[po(?)]- This relation follows with 



Vj~4^)(y) = /d2J'4^^(y,f)Vj~,po(f) 



„(2) 



(40) 



where CQ"'*(yi, . . . ,y„) is the n-particle direct correlation 
function in equilibrium, and 



from the more general form 

VjPo(?) + /3po(r)VjC/i(r) 

= Po(y1 fdv'c^^\f,f)Wf,po{?) 



(41) 



(42) 



of Eqs. (14) and (16) in reference [74]. Equations (38) 
and (|39|) lead to the equilibrium relation 



Spoil) 



(43) 



which we use instead of Eq. ( 37 ) as closure relation for 



Eq. ( 34 1 in the time-dependent (non-equilibrium) situa- 
tion, similar adiabatic approximation was used in the 
derivations of the DDFT equations for isotropic [TJ |3] 



and uniaxial |1] colloidal particles. The approximation 
results in the generalized DDFT equation 



-KAilt) 



dt 



(44) 



for anisotropic colloidal particles with the total equilib- 
rium Helmholtz free-energy functional 

F[pQ{f)] = J"id[po(yl] + -7'cxc[po(y)] + -7^cxt[/5o(y)] (45) 

that can be decomposed into the ideal rotator gas part 

m 



PF.Mf)] = /d2Jpo(y)(ln(A3po(y)) - 1) (46) 

with the thermal de Broglie wavelength A, the excess 
free-energy part -^oxc[po(y)], and the contribution [75] 



-Foxt[po(y1] = d^ poif)Ui{lt) 



(47) 



due to the external potential Ui{f,t). The DDFT equa- 
tion ( 44 ) describes the time evolution of the one-particle 



density for a system of similar anisotropic self-propelled 
colloidal particles that interact over a pair potential and 
is the main result of this paper. 



IV. SPECIAL CASES AND POSSIBLE 
APPLICATIONS 

There is no translational-rotational coupling in the uni- 
axial case, which means that I?"^^ (y^) and 'D^'^{f^) and 
therefore also I?™(ro) and I?f^(ro) vanish in this case. 
Furthermore, the one-particle density and the free-energy 
functional do not depend on the angle x for uniaxial par- 
ticles and the translational diffusion tensor can then be 
written as the matrix 'D'^'^{u) = D^^u^u + D±{l — u(E)u), 
which obviously only depends on the two independent 
short-time diffusion coefficients -Dy and D± for diffusion 
parallel and perpendicular to the orientation of the sym- 
metry axis u = (sin(6') cos(0), sin(0) sin(0), cos(0)) of the 
uniaxial particle, respectively, where 1 denotes the 3x3- 
dimensional unit matrix. Also the rotational diffusion 
matrix becomes quite simple for uniaxial particles. When 
we use T)^^ = DrI with the rotational short-time dif- 
fusion coefficient -Dr and the considerations above and 
neglect the self propulsion, we obtain the uniaxial DDFT 
equation from our more general DDFT equation ( 44 1 . 
From the uniaxial DDFT equation, one can in turn de- 
rive the DDFT equation for two spatial dimensions [7S] 
as well as the traditional DDFT equation for colloidal 
particles with spherical symmetry [1] |3] as special cases. 

The generalized dynamical density functional theory 
for passive and active biaxial particles as proposed in 



7 



Eq. (44) can be numerically solved for a plenty of dif- 



ferent problems. For passive particles, one can explore 
for example: i) the relaxation dynamics towards equilib- 
rium 0] , ii) the response of the system to time-dependent 
external potentials [77], iii) the growth of a thermody- 
namically stable phase into an unstable phase |78j . In- 
teresting effects for self-propelled particles include among 
others: i) the swarming and clustering behavior of biaxial 
particles in the bulk and in confinement [THl [75] , h) the 
combined impact of self-propulsion and external forcing 
[75] , iii) the effect of space- and time-dependent internal 
forcing |5D] . 



V. CONCLUSIONS AND OUTLOOK 

In conclusion, starting from the multi-body Smolu- 
chowski equation, we have derived dynamical density 
functional theory for self-propelled Brownian colloidal 
particles with arbitrary shape. This study was moti- 
vated by recent progress in synthesizing colloidal par- 
ticles with (almost) arbitrary shape. Our results consti- 
tute an important framework for further numerical explo- 
rations. This is in particular appealing as since recently 
an equilibrium density functional is known for arbitrar- 
ily shaped hard colloids [H] IHU [82] which can serve as an 
input for the dynamical density functional theory. An- 
other possibility to construct a density functional for bi- 
axial particles is a mean-field approximation for repulsive 
segment potentials @], which works for soft interactions 



[55] , or a perturbation theory jM] [SS] for anisotropic at- 
tractions around a spherical reference system. A large 
number of dynamical problems can then in principle be 
addressed including the dynamics [5BH55] and relaxation 
of nematic-like order in confined systems [U [SU] , nematic 
phases driven by external fields [77], nucleation kinetics 
of liquid crystalline phases [SUHSS] , and collective behav- 
ior of self-propelled particles [751 [75] . The results can 
be checked against Brownian dynamics computer simu- 
lations [94 1 195 ] . 

Possible extensions for the future are the inclusions of 
hydrodynamic interactions between the particles which 
are mediated by the solvent. Dynamical density func- 
tional theory of spherical particles was generalized for 
hydrodynamic interactions [55H5F. but this is not yet 
done for anisotropic particles. Another interesting exten- 
sion would be towards molecular dynamics which is the 
appropriate dynamics for molecular liquid crystals. But 
even for spheres it is much more complicated to formu- 
late a dynamical density functional theory for molecular 
dynamics [99Ul01j . Finally, the theory can readily be 
generalized towards binary mixtures [102j . 
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